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Abstract. Many conservative partial diflerential equations correspond to geodesic 
equations on groups of diffeomorphisms. Stability of their solutions can be studied 
I by examining sectional curvature of these groups: negative curvature in all sections 

implies exponential growth of perturbations and hence suggests instability, while 
positive curvature suggests stability. In the first part of the paper we survey what 
we currently know about the curvature-stability relation in this context and provide 
detailed calculations for several equations of continuum mechanics associated to 
Sobolev H° and energies. In the second part we prove that in most cases (with 
some notable exceptions) the sectional curvature assumes both signs. 

o\ 

Oh AMS Subject Classification (2000): 53C21, 58D05, 58D17. 

Keywords: Riemannian metrics, diffeomorphism groups, sectional curvature, stability, Euler- 
Arnold equations. 

g Contents 

1. Introduction [T] 

2. Preliminaries: metrics and geodesic equations [3] 
J> 2.1. The Euler- Arnold equations on Lie groups: examples [3] 

2.2. Further examples of infinite-dimensional geodesic equations [5] 

3. Global aspects of infinite-dimensional Riemannian geometry [S] 

3.1. Degenerate distance functions [S] 

3.2. Completeness and minimizing geodesies [S] 
O 3.3. Exponential map and extendability of geodesies [S] 
^ 4. Jacobi fields, curvature, and stability [IT] 

5. The sign of the curvature: previous results [IS] 

6. The sign of the curvature: the one dimensional case [IS] 
^ 7. The sign of the curvature: higher dimensions 120] 

^ 7.1. The F^-metric on Diff(T"): the EPDiff equation M 

7.2. The L^-metric on Diff(T"): the Burgers equation 122] 

7.3. The homogeneous i^^-metric on Diff^^exCT^) [23] 
Appendix A. Semi-direct products and curvature formulas [21] 
References [25] 



00 



1. Introduction 

The idea that stability of a dynamical system can be investigated using tools of 
Riemannian geometry goes back to Hadamard [22j who studied the free motion of 
a particle on a surface of constant negative curvature. A similar approach was con- 
sidered by Synge [65] . Perhaps the most influential example is due to Arnold [T] 
who showed that fluid motions can be viewed as geodesies in the infinite-dimensional 
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group of volume-preserving diffeomorphisms. This led him to examine curvature of 
the diffeomorphism group and derive a number of results on stability of ideal fluids. 
Roughly speaking, on a finite-dimensional manifold negative sectional curvature is 
related to instability, and positive sectional curvature is related to stability of the 
corresponding geodesic flow; the Ranch comparison theorem makes this comparison 
rigorous. Since the work of Arnold other partial differential equations have been 
interpreted as geodesic equations in infinite-dimensional spaces and, as with ideal 
fluids, calculating sectional curvatures in these cases has become a matter of broader 
interest. 

In this paper we examine the sign of the sectional curvature of certain metrics on in- 
finite dimensional manifolds (which are associated with several well-known equations 
of mathematical physics) and its relevance in the stability analysis of the associated 
initial value problems. 

Our main interest is in those equations that arise from right-invariant metrics on 
the group of (smooth) diffeomorphisms Diff(M), or its subgroup of volume-preserving 
diffeomorphisms (volumorphisms) Diff^(M), of a compact n-dimensional manifold M 
without boundary. Both spaces can be completed to Hilbert manifolds Diff*(M) and 
Diff^(M) modelled on Sobolev spaces of H" vector fields and divergence- free vector 
fields, respectively, with s > n/2 + 1. We assume that M has a Riemannian metric 
(•, •) with volume form ^. We equip the groups with right-invariant Sobolev metrics 
such that, on the tangent space at a diffeomorphism ry, we have 

(1.1) {{u o rj,v o r]))ri = j (^a{u,v) + b{6u^ ,6u^) + c{du^ ,du^)^ dfi, 

for any vector fields u, v on M. Here d is the exterior derivative, 5 = ziz*d* is its 
(formal) adjoint, a, b, c are non-negative constants and b and jj denote the standard 
"musical isomorphisms" of the metric corresponding to lowering and raising of indices. 
Formula (1.1) simplifies in dimensions n < 3 where we hav^ 



(1.2) {{uor],vo7]))n 



{auv + buxVx) dx n 

M 




+ b div u • div v + c curl u ■ curl vj dA n = 2, 
a{u, v) + 6 div u • div v + c(curln, curl v)) dV n = 3. 



Using the metric (1.1) one derives a number of PDE that are of interest in contin- 
uum mechanics and geometry. For example, in the one-dimensional case one obtains 
Burgers' equation, the Camassa-Holm equation, and the Hunter-Saxton equation as 
geodesic equations for appropriate choices of a, b, c. In higher dimensions one gets 
the EPDiff equation and the so-called template-matching equation. Furthermore, 
projection onto Diff^(M) yields the usual Euler equations of hydrodynamics and the 
Lagrangian-averaged Euler-a equation. We refer to the paper p9] and the references 
therein for more details. 

The goal of this paper is two-fold. In the first part we describe some aspects 
of Riemannian geometry of infinite-dimensional manifolds which are relevant to the 
analysis of partial differential equations of mathematical physics. In particular, we 
review the framework of the Euler-Arnold equations on Lie groups equipped with 
right-invariant metrics and explain the role played by sectional curvature in the study 
of (Lagrangian) stability. 

-'^Recall that curlu = *du^ and sgrad/ = —((5*/)" in two dimensions, while curlu = {*du^)'^ in 
three dimensions. 
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In the second part we present new results on the sign of the sectional curvature for 



the metric (1.1) on Diff(M) and Diff^(M) for different choices of the parameters a, b 
and c. In order to simplify calculations, we will take M to be either the circle or 
the flat torus T". We show that in most cases the sectional curvature assumes both 
signs. Two notable exceptions are: b = c = and n = 1, in which case the sectional 
curvature of Diff(S'^) turns out to be non- negative; and a = c = and n > 1, in 
which case the sectional curvature of the quotient space Diff(M)/Diff^(M) is strictly 
positive for any compact M. The latter case is studied in detail in the paper |29| . 
For the metric on Diff(5^) we obtain a simple curvature expression and explain 
how it can be viewed as the Gauss-Codazzi formula for an isometric embedding of 
the group of circle diffeomorphisms in a larger space. 



Besides the fact that right-invariant metrics (1.1) arise as Lagrangians of many 
PDE of continuum mechanics, our motivation to study them is also purely geometric: 
they are the natural metrics induced on orbits of pullback actions on spaces of tensor 
fields. For example, the canonical metric on the space of all Riemannian metrics 



on a compact manifold M induces a metric of the type ( 1.1 ) on Diff(M) (viewed as an 
orbit of any particular metric under the pullback action) whose geometry was studied 
in [29]. 

2. Preliminaries: metrics and geodesic equations 

In this section we review some well-known examples of infinite-dimensional Rie- 
mannian manifolds and their geodesic equations. Most interesting from our point of 
view are infinite-dimensional Lie groups, especially diffeomorphism groups equipped 



with right-invariant metrics. We describe examples of such groups in Section 2.1 
In Section 12.21 we discuss other situations of interest when the manifolds are not Lie 
groups or the metrics are not right- invariant. 

2.1. The Euler- Arnold equations on Lie groups: examples. We begin by de- 
scribing the general Lie-theoretic setup of Arnold [Ij. 

Consider an infinite-dimensional Lie group G equipped with a smooth right-invariant 
(weak) Riemannian metric determined by an inner product ((•,•)) defined on the tan- 
gent space TgG at the identity element. A geodesic in the group starting from e in 
the direction uq can be obtained from the solution u{t) of the Cauchy problem for 
the associated Euler- Arnold equation on T^G, namely 

dzL 

(2.1) — = -ad>, u{0) = no 

where u — )• ad*u is the adjoint of the linear operator u — )• adyU = —[v, u] with respect 
to the inner product on T^G, that i^ 

(2.2) {{adlu,w)) = -{{u,[v,w])), u,v,w£TeG. 
The geodesic is now obtained by solving the flow equation 

(2.3) ^ = DR^(^t)u{t), m = e 
where ^ i— RriiO denotes the right-translation in G by ry. 



^Here [•,•] denotes the commutator on T^G induced by the Lie bracket of right-invariant vector 
fields on G, i.e. [u, v] = [X, Y]^ where X, Y are the right-invariant vector fields determined by X^ = it 
and Ye = v. 
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In the case of ideal hydrodynamics G is the group of volumorphisms (volume- 
preserving diffeomorphisms) of a manifold M 

Diff^(M) = {r/ G Diff(M) | = 

with the right-invariant metric given at the identity by the inner product, i.e. by 



setting 6 = c = 0in(l.l). The resulting Euler- Arnold equations (2.1 ) are the familiar 



Euler equations of incompressible fluids in M 



(2.4) 



ut + V„u = — Vp, divtt = 0, 



frequently written in the form 

(2.5) ut + CuUJ = 0, 

where uj = curlu is the vorticity]^ 

If we set a = 1 and c = with 6 = in ( 1 . 1) on the volumorphism group then the 



corresponding Euler- Arnold equation (2.1) is called the Lagrangian-averaged Euler 
equation [Ml EH- 
curlu 



UJ 



It is more complicated than (2.4), but its "vorticity" given by 
a^Acurln satisfies the same equation as (2.5). The analysis of this 



equation presents similar difficulties as (2.4) or (2.5); see e.g. 



Remark 2.1. The volumorphism group is one of the three "classical" diffeomorphism 
groups. Another is the symplectomorphism group Diff(j(M), consisting of diffeomor- 
phisms preserving a symplectic form a; of an even-dimensional manifold M, and the 
third is the contactomorphism group DiffQ,(M) consisting of diffeomorphisms r] such 
that r]*a = Fa, where a is a contact form and -F is a nowhere-zero function on an 
odd-dimensional manifold M. Geodesic equations of the right-invariant metric on 
these groups have been studied in [T31 dSl [TB] . In two dimensions the geodesic equa- 



tion on the symplectomorphism group reduces to the Euler equation (2.4), and in one 



dimension the geodesic equation for contactomorphisms reduces to the Camassa-Holm 
equation. Simpler equations arise when the metric is restricted to the subgroup 
of Hamiltonian diffeomorphisms or the subgroup of strict contactomorphisms; see 
Smolentsev [64J for a review of their properties. 

It is worth pointing out that the subgroup of Hamiltonian diffeomorphisms carries 
a bi-invariant metric given at the identity by 



(2.6) 



((sgrad /, sgrad g)^ 



M 



where / and g are assumed to have mean zero. For such metrics we have {{u, ad^w)) + 
{{advU, w)) = whenever u, v and w are in the Lie algebra. It follows that ad* it = 0, 
and hence the Euler equation (2.1) reduces to du/dt = 0. Geodesies are easy to find 
as they are one-parameter subgroups: simply fix a velocity field and compute the 
flow. 

The group of circle diffeomorphisms Diff(5^) has been a rich source of examples. 
In this case we can set c = in the formula ( |1.1[ ) since in one dimension dv!' = for 
any vector field u. Two much-studied Euler- Arnold equations that arise here are the 
periodic (inviscid) Burgers equation 



(2.7) 



ut + 3uux = 



In two dimensions ui 
a vector field and £„cj = [u, lo] 



curl M is a function and 



{u, Voj); in three dimensions u: = curltt is 
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associated with the inner product (with 6 = 0) and the periodic Camassa-Hohn 
equation 

(2.8) Ut — Utxx + "iuUx — 2UxUxx — UUxxx = 

obtained from the product (with o = 6 = 1) [3l 138]. 

Interesting examples also arise on the Bott-Virasoro group Vir(S'^), the universal 
central extension of Diff(5^), with a group law defined by 

(2.9) (r?,a)o(e,/3) = (^r, o + (3 + ^ j^^logdxiv o dlogdx^^ 

where ??, ^ S Diff(S'^) and a,/3 G M. The right-invariant metric on the Bott-Virasoro 
group given at the identity by the inner product 

{{{u,a),{v,b)))]^2 = / uvdx + ab, 



yields as its Euler- Arnold equation (2.1) the Korteweg-de Vries equation [51 152| 
(2.10) Ut + 3uux + auxxx = 0, a = const. 

If we use the inner product instead 

{{{u,a),{v,b)))Hi = / {uv + UxVx) dx + ah, 



then the corresponding Euler- Arnold equation is the Camassa-Holm equation with 
drift \48i 

Ut — Utxx + 'iuUx — 2UxUxx + {k — u)Uxxx = 0, K = COUSt. 



Equations (2.7) and (2.8) have higher-dimensional analogues. The n-dimensional 
version of the Burgers equation (2.7) arising from a right-invariant metric on 



Diff(M) is the so-called template-matching equation [5T] 

Ut + VuU + {Vu)\u) + (divii)M = 



while the n-dimensional version of the Camassa-Holm equation (2.8) on Diff(M) is 
the EPDiff equation [S^. 

Other examples of geodesic equations for right-invariant metrics on infinite-dimensional 
groups include 

• the equation of passive scalar motion (on the semidirect product of Diff^(M) 
with C~(M)); 

• the equation of 3D magnetohydrodynamics (on the semi-direct product of 
Diff^(M) with divergence- free vector fields) 

ut + VuU = -Vp + VbB, Bt + [u,B]=0, divii = divB = 0; 

• the fi-CH equation, corresponding to an H^-type metric on Diff(S'^), see Re- 



mark 6.4 



• the quasigeostrophic equation in /3-plane approximation on T'^ on the central 
extension of the group of Hamiltonian diffeomorphisms 

Ut + uj} = -fidxip, uj = A^; 

as well as the Boussinesq approximation to stratified fluids, equations for charged 
fluids and fluids in Yang- Mills fields, the Landau-Lifschitz equations as well as various 
2-component generalizations of the one-dimensional equations mentioned above; see 

e.g., mm- 

2.2. Further examples of infinite-dimensional geodesic equations. 
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2.2.1. Spaces of curves. Let VtM be the (free) loop space over a compact Riemannian 
manifold M whose points are smooth maps from to Mj^ The tangent space to 
VlM at a point 7 consists of vector fields in M along 7, i.e., maps s — )• V{s) G T^(g-)M. 

Two metrics on VtM have been of particular interest. The first is the weak Rie- 
mannian !?■ metric given at 7 G Q.M by 



(2.11) {{U,V))l2,^= {U{s),V{s))^^,)ds. 



Its geodesies correspond to geodesies on the underlying manifold: if r/(0) = 7 and 
77(0) = V, then 7]{t){s) = exp^(,)(iy(s)). 
The other is the Sobolev metric 



ds, 



(2.12) HU. V))„., = I [msl V(S)U, + ^)_^^^.^ 

where D/ds denotes the covariant derivative along 7 in M. The metric (2.12) is in 
fact the more natural of the two and, in particular, turns the set fiM^^ consisting of 
all loops in M into a complete Hilbert Riemannian manifold. 

The set of simple closed curves in can also be regarded as a version of the 
loop space. Such curves may be viewed as boundaries of planar "shapes" and finding 
a suitable notion of distance between "shapes" has been of interest in applications 
to pattern theory. For this purpose, however, parameterizations are irrelevant; thus 
it is useful to pass to the quotient by the diffeomorphism group of . Geodesies 



on the quotient of the metrics induced by (2.11) and (2.12) were studied by Michor 
and Mumford [l3] . Another approach is to consider the subspace consisting of those 
curves parameterized by arc length (or its multiple) with the induced or metric. 
The geodesies on this subspace are solutions of a wave-like equation 

(2.13) rjtt = dsiaijs) , cr^s - h^spcr = -|r/stp, |r/s| = 1, 

which describes an inextensible string (the asymptotic limit of a string with a very 
strong tension) ; see [58j for a geometric discussion of this equation. 

2.2.2. Homogeneous spaces. Degenerate metrics on diffeomorphism groups also lead 
to natural geometries on their quotient spaces. For example, if it is a vector field on 
and rj G Diff(S'^), then setting 



{{u ori,vor]))^i = / Uj:Vx dx 



one obtains an invariant degenerate Sobolev metric which is a limiting case of 
the metric (1.2) when a — )• or 6 — )• 00. It becomes a weak Riemannian metric 



on the homogeneous space Diff(5^)/S'^. The corresponding Euler- Arnold equation is 
the Hunter-Saxton equation [S^ 

(2.14) Utxx + '^UxUxx + UUxxx = 0. 

Passing to the quotient space is geometrically appealing since the manifold turns out 



to be isometric to a subset of the round Hilbert sphere 

More generally, the same construction applies on any Riemannian manifold M 
using the right-invariant degenerate metric 

{{u o ?y, v o T])) £fi = / div u ■ div v dfi, 



^If Af is a Lie group then ilAf becomes a loop group under pointwise multiplication. 
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on the quotient space Diff(M)/Diff^(M). Its Euler- Arnold equation is a higher- 



dimensional analogue of (2.14) given by 



(2.15) Vdivut + divnVdiv'u, + V(u, Vdivu) = 0, 

and one can establish a similar isometry with the round sphere in a Hilbert space. The 
induced Riemannian distance turns out to be a spherical analogue of the Hellinger 
metric in probability theory; see [29j for details. 

2.2.3. Spaces of maps and non-invariant metrics. More generally, given a Riemannian 
manifold M and a compact manifold N with a volume form v (and possibly with 
boundary) consider the space C°°{N,M) of smooth maps from N into M. On each 
tangent space at / G C°°{N,M) we can define an metric by 

(2.16) iU,V))f= [ {U{x),V{x))f^^)du{x), 

J M 



and, as in the case of the loop space and (2.11 ), show that its geodesies come directly 
from geodesies on M. 

The group of smooth diffeomorphisms Diff(M) is an open subset of the Frechet 



manifold C°°(M, M) so that we can likewise put the metric (2.16) on it. Note that this 
metric is not right-invariant; nevertheless the corresponding geodesic equation can be 
rewritten on the tangent space to the identity where it becomes the multidimensional 
inviscid Burgers (or pressureless compressible Euler) equation 

(2.17) _+V„u = 0. 

Since geodesies in M starting from two nearby points will inevitably cross (at which 
time the geodesic in Diff(M) must exit the diffeomorphism group) solutions of the 
pressureless Euler equation solutions in general will blow up in finite time. Physically, 
this corresponds to the emergence of a shock wave leading to collisions of the fluid 
(or gas) particles. Nonetheless, the geodesic remains in C°°(M, M) for all time. 

Remark 2.2. The equations of incompressible fluids with boundary can be viewed 
formally as geodesic equations on the space £^{Q.,W^) of volume-preserving embed- 
dings of the closure of an open subset of into W^. These equations were studied 
geometrically by Ebin pTS] and shown to be identical to the standard equations of 
incompressible fluid mechanics except for the fact that the boundary condition for 
the pressure is a Dirichlet rather than a Neumann condition. 

2.2.4. Spaces of metrics. Another geometrically interesting space is the space Met(M) 
of all Riemannian metrics on a compact manifold M. For any metric g on M and any 
77 € Diff(M) we define the pullback metric r/*g. In jllj Ebin studied the Riemannian 
metric on Met(M) which is right-invariant under the pullback action. Given g G 
Met(M) and tangent vectors A and B (smooth tensor fields of symmetric bilinear 
forms) the metric is defined by 



{{A,B))g= [ TTg{AB)dl^g, 

Jm 



where Tr^ is the trace with respect to g and fig is the Riemannian volume form. 
The curvature and geodesies of this metric were computed explicitly by Freed and 
Groisser [19] : sectional curvature is non-positive and geodesies generally exist only for 
finite time (until the metric becomes degenerate). The diffeomorphism group embeds 
in Met(M) as an orbit of a generic g (i.e. with no non-trivial isometrics) and if g is 
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Einstein then the induced metric on Diff(M) is a special case of (1.1). We refer to 
Clarke [8j for recent results on the distance and diameter of this space. 

Similarly, one can endow the space of all volume forms Vol(M) on M with a natural 
right-invariant metric given for n-form fields a and (3 tangent to ^ by 

{{a,(3))u = f --dv. 

Although this metric is flat, it too is not geodesically complete in general. Orbits 
of the diffeomorphism group in Vol(M) are the homogeneous spaces of densities 
Diff(M)/Diff^(M) of constant positive curvature; we refer to [29j. 

The map g fig where fig is the Riemannian volume of 5 is a submersion which 
becomes a Riemannian submersion after suitable rescaling of the metric on Vol(M). 
Its fibers are the spaces of metrics Met^(M) with the same volume form /x. These 
fibers are globally symmetric (with negative curvature and indefinitely extendable 
geodesies) in the induced metric from Met(M). The natural action on Met^(M) is 
pullback by volumorphisms, and if g has no nontrivial isometries, then the orbits of 
Diff^(M) are embedded submanifolds with right-invariant metrics. 



3. Global aspects of infinite-dimensional Riemannian geometry 

The obstacles that arise in the study of global Riemannian geometry of infinite 
dimensions manifolds are well known. They are mostly caused by the lack of local 
compactness or the fact that the topology generated by the metric may be weaker 
than the manifold topology. As a result some of the finite-dimensional techniques are 
not available or are of limited use. For example, the Riemannian exponential map 
may not be defined on the whole tangent bundle or even be smooth, conjugate points 
may cluster along finite geodesic segments or have infinite multiplicity, etc. In this 
section we illustrate some of these situations with a few familiar examples. 



3.1. Degenerate distance functions. The distance between two points in a weak 
Riemannian Hilbert manifold can be defined as in finite dimensions, i.e., as the in- 
fimum of lengths of piecewise smooth curves joining them. It is easy to prove that 
it satisfies all the axioms of a metric space except for nondegeneracy which typically 
requires some additional assumptions. 

As an example, consider the right-invariant metric on Diff(M) defined by 



{{uo r],v o ' 



{u, v) dfi, 



M 



which corresponds to the case 6 = c = in ( 1 . 1 ) . In |42| it is shown that the geodesic 



distance is identically zero, (i.e., between any two diffeomorphisms there are curves 
of arbitrarily short length). This is essentially related to the lack of control over the 
Jacobian. The same phenomenon also occurs for the right-invariant metric on the 
Bott-Virasoro group [5j, for the metric on the "shape space" of curves modulo 
reparameterizations [l3], and for the bi-invariant Hofer-type metric (2.6) on the 
Hamiltonian diffeomorphisms |17j. 

On the other hand, we obtain nondegenerate Riemannian distances for the 
metric on Diff^(M) [H], the right-invariant metric on Diff(M) corresponding to c = 
in (1.1) [50], the space of maps C°°{N^M) in the metric (2.16), the space of arc- 



length parameterized curves [58j, and the space Diff(M)/Diff^(Af) of densities p9] . 
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3.2. Completeness and minimizing geodesies. Even if the geodesic distance is 
nondegenerate, thus providing a genuine metric space structure on a space of maps, 
this metric space may not be complete. For example, the completion of {N, M) 
in the non-invariant metric (2.16) consists of measurable maps from to M 



which may not even be continuous. The same phenomenon occurs for the group of 
volumorphisms Diff^(M) with the metric if M is a three-dimensional manifold — 
in this case the completion in the Riemannian distance is the space of all measure- 
preserving maps; see [62] 

In finite dimensions, completeness of a Riemannian manifold M as a metric space is 
equivalent to geodesic completeness, i.e., extendability of geodesies for all time, which 
in turn implies that any two points in M can be joined by a minimal geodesic. The 
proof of this result, the Hopf-Rinow theorem, relies crucially on local compactness, 
and the result is no longer true in infinite dimensions as observed by Grossman [21j 
and Atkin |3] . The former constructed an infinite-dimensional ellipsoid in the space 
of square-summable sequences with points which cannot be connected by a minimal 
geodesic in the induced metric from and the latter modified this construction to 
get points that cannot be joined by any geodesic at all. Interestingly, Ekeland |16j 
showed that on a complete Riemannian Hilbert manifold the set of points attainable 
from a given one with a minimizing geodesic contains a dense Gs (i.e., a countable 
intersection of open sets). 

One situation in which everything works nicely is the completion of the space 
of smooth loops QM^^ . Unlike many other examples discussed here, it is a genuine 
(strong) Riemannian Hilbert manifold in the topology generated by the distance 



function of (2.12), and by the result of Eliasson [18] any two of its points can be 



joined by a minimizing geodesic. 

3.3. Exponential map and extendability of geodesies. If the Cauchy problem 
for the geodesic equation on a (possibly weak) Riemannian Hilbert manifold M is 
locally well-posed, then the exponential map of M can be defined as in finite dimen- 
sions. Using the scaling properties of geodesies we set 

(3.1) expp-.U cTpM-^ M, expp(i;) = 7(1), 

where 7(i) is the unique geodesic from 7(0) = p with initial velocity 7(0) = v in some 
open neighbourhood U of zero in the tangent space at p. 

In general local (in time) well-posedness refers to constructing a unique solution for 
a given initial data on a short time interval which depends at least continuously on 
the data. However, from the point of view of difi^erential geometry, it is desirable if the 
dependence is at least smooth. Indeed, in this case applying the inverse function 
theorem for Banach manifolds, it is possible to deduce that (as in finite dimensions) 
the exponential map is a local diffeomorphism; this implies in particular nondegen- 



eracy of the geodesic distance as in Section 3.1, Furthermore, other geometric tools 
such as Jacobi fields and curvature can be introduced to study rigorously stability as 
the problem of geodesic deviation (we shall elaborate on this in Section [4]). 

The exponential maps defined on suitable Sobolev completions in the examples dis- 
cussed so far are either at least smooth or else are continuous (even differentiable) 
but not C^. The former include 

• the metric on the volumorphism group Diff^(M) whose geodesies corre- 
spond to the Euler equations of ideal hydrodynamics; see [2], 



^If M is two-dimensional, the completion of T>iS^{A4) in the metric is unknown. 
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• the metric on Diff^(M) corresponding to the Lagrangian-averaged Euler 
equation; see [6T] . 

• the metric on Diff''(5^) corresponding to the Camassa-Holm equation in 
[2] and its generahzation to Diff^(M) and the EPDiff equation in |50j . 

• the metric on the free loop space QM"^ in |Tfj, 

• the homogeneous metric on Diff'^(S'^)/5'^ corresponding to the Hunter- 
Saxton equation in |30] and its generahzation (2.15) on DifF^(M)/Diff* (M); 
see [29], 

• the right-invariant metrics on Met'*(M), Met^(M) and VoP(M) in [TUfTO]. 

• the noninvariant metric on H^{N,M) or Diff'^(M) whose geodesies are 
described by pointwise geodesies on M; see [11]. 

The metrics for which dependence fails include 

• the metric on the Virasoro group whose geodesic equation corresponds to 
the Korteweg-de Vries equation; see [lOj, 

• the right-invariant metric on Diff'^(S'^) which yields the (right-invariant) 
Burgers equation [9] or its higher-dimensional generalization and the template- 
matching equation |42j . 

• the metric on unit-parametrized curves in the plane yielding the whip equa- 
tion (2.13) in t58j, or on the equivalence classes of curves under reparametriza- 



tions in |43j . 

Later on we will describe examples where the exponential map fails to be as a 
result of accumulation of conjugate points at i = (as in [9], [10], [l6], and |58|). 

It is well-known that the group of volumorphisms Diff^(M) of a two-dimensional 
manifold M equipped with the right- invariant metric is geodesically complete, i.e., 
its geodesies which correspond to solutions of the 2D incompressible Euler equations 
are defined globally in time when s > 2|^ This result is due to Wolibner [68] with sub- 
sequent contributions by Yudovich ^69j and Kato [27 \ and follows from conservation 
of vorticity, although the argument is not routine. In three dimensions the problem 
is open and challenging. 

One might expect the right-invariant metric on the volumorphism group to be 
somewhat better behaved but as of now we have the same result: global existence in 
two dimensions is relatively easy ^61j but in three dimensions is unknown |25j . On the 
other hand, it is known that smooth solutions of the one-dimensional Camassa-Holm 
equation (2.8) break down for certain initial data |40ll41j . All solutions of the periodic 
Hunter-Saxton equation (2.14) as well as its higher-dimensional generalization (2.15) 
are also known to blow up in finite time. 

For the non-invariant metric on DifF*(M) (whose geodesies are given by point- 
wise geodesies on M) global existence clearly fails: for a typical initial velocity field 
two geodesies will eventually cross (which corresponds to a "shock"). On the other 
hand in H^(M, M) geodesies exists for all time since such maps need not be injective. 

Geodesies in the space of metrics Met*(M) and volume forms VoF(M) typically 
become degenerate in finite time, see [8] and [19], while geodesies in Met* (M) persist 
for all time [11] . Both the and metrics on the homogeneous space of equivalence 
classes of curves in the plane admit geodesies that degenerate to points [43]. On the 
space of unit-speed curves, the whip equation would be expected to blow up in finite 
time physically [58], but this is not yet proved. 



We do not discuss weak solutions (for which much of the geometry seems to break down) of the 
PDE mentioned above. 
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4. Jacob: fields, curvature, and stability 

As mentioned in the Introduction, one application of Riemannian techniques in 
the study of equations of fluid dynamics has been to the problem of (Lagrangian) 
stability using the equation of geodesic deviation (the Jacobi equation) which involves 
the curvature tensor. In this section we describe this approach for a general infinite- 
dimensional manifold equipped with a possibly weak Riemannian metric but whose 
exponential map is assumed to be at least C^. We will discuss Jacobi fields (as 
infinitesimal perturbations) and the role played by sectional curvature and its sign. 
Various results for specific examples mentioned in the previous sections will be the 
subject of Section [5} 

Let be a (possibly weak) Riemannian Hilbert manifold whose geodesic equation 
is written in the form 

Ddj _ 
dt~dt ~ ' 

where ^ is the tangent vector field and ^ is the covariant derivative along the curve 
7(t) in A^. If j{s,t) is a family of geodesies with 7(0, t) = 7(t) then the formula 

(4.1) j(t) = |l(o,t) 

gives a Jacobi field J{t) along 7, i.e. a solution of the Jacobi equation 

(4.2) ?^ + <^>?)? = 
^ ' dt^ \ dtJ dt 

obtained by differentiating the geodesic equation in s and evaluating at s = 0. As in 
finite dimensions the Riemann curvature tensor R of M arises here due to the fact 
that covariant derivatives do not commute in general. Furthermore, the basic result 
of Cartan applies as well so that for any v and w G TpM we have 

(4.3) {De^Y>p)tv{tw) = J{t) 



where J{t) is the Jacobi field along 7(t) = expp(ti;) solving (4.2 ) with initial conditions 
J(0) = and J'(0) = w. 

Recall that the sectional curvature in the direction of the 2-plane spanned by the 
vectors Xp and Yp S TpAi is given by the formula 



(4.4) K{p) 



{R{Xp, Yp)Yp, Xp) 



Example 4.1. On a two-dimensional Riemannian manifold the Jacobi equation can 
be reduced to a single ODE for a function j(t) representing the component of J 
orthogonal to 7, which takes the form 

(4.5) f01+Ki^{t))j{t) = O, 

where K is the sectional curvature at a point 7(t). In the special case where K is 
constant the solution of (4.5) with j(0) = is 



j(i)=/(0) 



4^ sin VKt K >0 

t K = . 

sinh ^AKlt K <0 



/\K\ 

This simple special case is the source of much of our intuition about curvature and 
stability. Suppose that we know precisely the initial position of a particle traveling 
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along a geodesic and its initial velocity only approximately. If > then all Ja- 
cobi fields are bounded uniformly in time, and thus geodesies starting with nearby 
initial velocities will remain nearby for all time. If K < then the Jacobi fields grow 
exponentially in time, so that small errors are magnified and the motion is unpre- 
dictable, li K = then the growth is polynomial. On higher dimensional manifolds 
with variable curvature things become more subtle. 

Recall that singular values of the Riemannian exponential map are called conjugate 
points. More precisely, two points p and q along a geodesic in Ai are conjugate if 



D expp viewed as a linear operator from Tp^A to TqM. given by (4.3 ) either fails to be 
injective (in which case the points are called mono-conjugate) or it fails to be surjective 
(the points are called epi- conjugate); see Grossman [21j. In finite dimensions both 
types coincide. 

Next, we state the Ranch comparison theorem for weak Riemannian metrics follow- 
ing Biliotti This result relates growth of Jacobi fields to bounds on the sectional 
curvature and is a far-reaching generalization of Sturm's comparison theorem on os- 
cillation of solutions of second order ODE. 

Theorem 4.2. Let Ai and Ai be (possibly infinite- dimensional) weak Riemannian 
manifolds modeled on Hilhert spaces E and E, with E isometric to a closed subspace of 
E. Assume that A4 and A4 have smooth Levi-Civita connections (and hence smooth 
exponential maps) with sectional curvatures K and K . Let 7 and 7 be two geodesies 
of equal length and suppose that for every X G T^(()7W and X G Ty^f-^M 

K{X,j'{t))<K{X,^'{t)). 

Let J and J be the Jacobi fields along 7 and 7 such that 

• J(0) = and J(0) = 0, 

• J'(0) is orthogonal to 7'(0) and J'(0) is orthogonal to 7'(0), and 
^ . ||J'(0)|| = ||J'(0)||. 

If J{t) is nowhere zero in the interval (0, a] and if 7 has at most a finite number of 
points which are epi-conjugate but not mono-conjugate in (0,a], then 

(4.6) ||J(t)|| > ||J(t)|| forallte[0,a]. 

It often happens that such pathological points which are epi-conjugate but not 
mono-conjugate can fill out a whole interval (this is the case for the volumorphism 
group of a three-dimensional manifold [57]) so that the criterion above may only be 



useful if the exponential map is Fredholm; see Remark 4.6 below 
Remark 4.3. Theorem 



4.2 



implies that if K{X, 7) > —k for some positive constant 
k (i.e., take M to be a constant negative curvature space) and if 7 is free of conju- 
gate points, then any Jacobi field along 7 in M satisfies ||J(t)|| < \\J' {0)\\k-^/^ sinh kt, 
which gives essentially the maximum Lyapunov exponent for the system. In the oppo- 
site direction, if i^:(X,7) < (i.e., take M to be a flat space) then ||J(t)|| > || J'(0)||i, 
which can be interpreted as a weak instability with perturbations growing at least lin- 
early in time. In general, however, one should be cautious when drawing conclusions 
based on the Ranch theorem and one's finite-dimensional intuition: positive curva- 
ture does not imply stability, while negative curvature does not necessarily imply 
exponential instability, as we discuss below. 

For a general Riemannian manifold without additional structure one does not ex- 
pect more precise results on the relation between curvature and stability. However, 
most of our examples have a group structure under which the Riemannian metric is 
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right-invariant, which can be used to get additional information. To this end it will 



be convenient to decouple the Jacobi equation (4.2) into two first-order equations. 



Namely, let be a group G with a right-invariant (weak) Riemannian metric. As 



in Section 2.1 its geodesies j{t) can be described by a pair of equations consisting of 



the Euler- Arnold equation (2.1) and the flow equation (2.3) 

du 



(4.7) 



dj 
~di 



DRyiu), 



dt 



+ ad,,ti = 



defined in G and T^G respectively. Let j{s, t) be a family of such geodesies with 
7(0, t) = 7(t) and with Eulerian velocity u = DR^-i^. Setting y = DR^-iJ (where 



J is the Jacobi field along 7 as in (4.1)) and z = §^|s=o a-nd differentiating both 
equations in (4.7) with respect to s at s = we obtain a splitting of the Jacobi 



equation (4.2) into 

(4.8) 

(4.9) 



dy 

dt 
dz 

dt 



aduV 



+ ad*2; + ad*u = 0. 



The linearized Euler equation (4.9) can be used to define a notion of stability: a 
solution u of the Euler equation is (linearly) stable if every perturbation z is bounded 
uniformly in time. Using (4.8) one can then relate this notion to stability of La- 
grangian trajectories (as in |53j for the volumorphism group) and draw sharper con- 
clusions about the behaviour of geodesies than is generally possible using only Ranch's 
theorem. The next two examples illustrate the subtleties. 

Example 4.4 (Rigid body motion). Let G = SO{3). Its Lie algebra so(3) = T^G is 

spanned by the vectors ei, 62 and 63 satisfying [ei, 62] = 63, [e2, 63] = ei, [63, ei] = 62. 
This is the group of antisymmetric matrices represented as 



x^ei + x^e2 + x^e3 



The shape of the rigid body determines a left-invarianl|^ Riemannian metric by the 
conditions (ei,ei) = Ai, (62,62) = A2, and (63,63) = A3 for some positive numbers 






-x^ 


— X 







— X 




x^ 






Ai,A2,A3. The left- invariant analogue of (|4.7|) reads as follows 

ad*M 



dj du 

with the Euler- Arnold equation given explicitly by 

du^ A2 — A3 2 3 du"^ A3 — A 



A 



dt 



A2 



du^ 
"dt 



Ai - A 



A3 



2 1 2 
u u . 



Consider one steady solution of this equation given by u^ = 



with 



1, 



supposing that < Ai < A2 < A3. The linearized Euler equation (4.9) takes the form 
dz^ 



dt 



Ai 



df_ 
dt 



0, 



dz^ 
~dt 



Ai - A 



A3 



'z^ 



of our equations so far which have been stated for right-invariant metrics apply to left- 
invariant metrics after a possible change of the sign. 
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and if Ai < A2 < A3 it admits exponentially growing solutions]^ The linearized flow 



equation (4.8) takes the form 



dy^ , 3 1 dy'^ 2 dy^ 1 3 

— \-y = z , —r- = z , — y = z , 

at at at 

so that y{t) grows exponentially if z[t) does. On the other hand, sectional curvature 
in all directions containing 62 can be made positive by a suitable choice of A's. Using 
e.g., the formulas of Milnor [41] we compute for x = x^ei + x^e2 that 

in( N \ (A2 - Ai)^ - A| + 2A3(A2 - A3 + Ai) 12 



4A 



3 



(A3 — A2)^ ~ -^1 + 2Ai(Ai — A3 — A2) , 3.2 
+ 4Ai ' ' 

which is positive-definite if e.g., Ai = |, A2 = 1, A3 = |. Hence, we have positive 
curvature along the geodesic but exponentially growing Jacobi fields. This happens 
because the Ranch comparison theorem bounds Jacobi fields only up to the first 
conjugate point; beyond that point all bets are off. 

Example 4.5 {Couette flow). Consider M = Diff^([0, 1] x S^) with the metric. 
The "plane-parallel Couette flow" u{x, y) = x is a. steady solution of the Euler 
equation. This solution is known to be Eulerian stable [21 [53] even though the sec- 
tional curvature is non-positive in all sections and typically negative |45j . A closer 
inspection reveals that the growth of all Jacobi fields is precisely linear [5l]. Hence, 
we do no better than what Ranch's theorem says: negative curvature need not imply 
exponential instability. 

It is worth pointing out that one can relate Eulerian stability to Lagrangian 
stretching, i.e., to the growth of Ad^(4) in the operator norm. Using the formula 
d{Ad.f^y) / dt = Adr]{dy / dt — aduy) and defining Y and Z hy y = Ad^y and z = AdjjZ, 
equation (|4.8l) can be rewritten as 



while equation (4.9) becomes 



(4.11) ^(Ad;Ad^Z) + ad>o = 

after incorporating conservation of vorticity Ad*n = uq; see |50]. Observe that the 
operator Z — )• Ad*^j^Ad^(()Z is selfadjoint and positive-definite while Z 1— )• ad^uo is 



anti-selfadjoint with constant coefficients. This makes (4.11) somewhat simpler to 



analyze than (4.9) (even if u is independent of time) because z 1— )■ ad*z -|- ad*M is not 



selfadjoint; see [63] . 

Remark 4.6 {Fredholm exponential maps and conjugate points). If the exponential 
map of a weak Riemannian manifold A4 is known to be smooth, one can ask about the 
distribution and nature of its singular values (conjugate points): can mono-conjugate 
and epi-conjugate points coincide, have finite multiplicity, or be discretely distributed 
along finite geodesic segments? 



This corresponds to the well-known fact that rotation of a rigid body about its largest and 
smallest axes is stable but rotation about the middle axis is unstable. 
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These questions turn out to have positive answers if the derivative of the exponen- 
tial map is a Fredholm operator between the tangent spaces of A4 with index zeroj^ 
This was first established for the free loop space QM"^ with the Sobolev metric 



(2.12) in |47j . In this case the proof of Fredholmness is based on the fact that the 



curvature operator R in (4.2) is compact. In general one does not expect compact- 
ness. However, in the special case when is a group G with a right-invariant metric 
one can analyze the derivative of the exponential map using the pair of equations 



(4.10)-(4.11 ) to conclude that it is a sum of two terms, the first determined by the 
positive-definite operator Z — )■ Ad* Ad,,Z and the second a composition of a bounded 
map with Y — t- adytiQ. The former is invertible. Thus, if the latter is compact then 
the exponential map will be Fredholm of index zero. This strategy works on the 
volumorphism group Diff^(M) with the metric if n = 2 (but not if n > 3) as well 
as for right-invariant Sobolev metrics of sufficiently high order on any diffeomorphism 
group; see [50] . 

Fredholmness fails for the exponential map of the ellipsoid in with the in- 
duced metric (there are sequences of mono-conjugate points accumulating at an epi- 
conjugate point, or a mono-conjugate point of infinite order [21]). It also fails for the 
free loop space with the metric |47| and for the volumorphism group Diff^(M) of 
a three-dimensional manifold in the metric |57j . In the latter case, we can have 
mono-conjugate points that are dense in an interval and epi-conjugate points that fill 
up an interval. 

On the other hand, smoothness of the exponential map implies (by the inverse 
function theorem) that any sufficiently short geodesic segment is free of conjugate 
points and hence is locally minimizing. In particular, if for some tn \j the points 
j{tn) along a given geodesic are mono-conjugate to 7(0) then the exponential is not 
C^. This method was used to prove that the exponential maps associated to the 
KdV equation [46J, the right-invariant Burgers equation [9J, and the whip equation 



(2.13) [58j cannot have exponential maps. 

5. The sign of the curvature: previous results 

In the remainder of the paper we will focus on the sign of the sectional curvature 
in the examples described above. It turns out that with few exceptions sectional 
curvature can be positive or negative depending on the two-dimensional direction. 
This section contains a survey of known results and techniques. New results will be 
presented in Sections [6] and [7} 

The simplest curvature formula arises on Diff(M) C C°°(M, M) equipped with the 
metric KTol). If U = uo r] and V = v o r] are two vector fields on Diff(M) where 
n, f G TeDiff(M) then the covariant derivative of the metric is computed in terms 
of the covariant derivative on M as {'VuV)'^ = (Vuf)ory. Therefore, the curvature 
of Diff(M) is completely determined by the Riemannian curvature of M and can be 
computed directly from the definition 

{{R{u o ri,v o ri)v o ri,u o 1])) 1^2 = / {R{u,v)v,u) o 7] dfi. 

Jm 

A few simple observations can be made based on this formula. If the vector fields u 
and V have unit norms and are chosen to have disjoint supports then the integral 
on the right-hand side will be zero. Consequently, the curvature of Diff (M) cannot 
be strictly positive or strictly negative even if M has constant (non-zero) curvature. 

'^In this case the exponential map is said to be a nonhnear Fredholm map of index zero (provided 
that A4 is connected). 
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Furthermore, it is also clear that it can be non- negative (non-positive) if and only if 
the sectional curvature of M is non- negative (non-positive). 

One can similarly obtain relatively simple expressions for the curvature of the 
volumorphism group Diff^(M) (following |l5]) or the free loop space ilM as well as 
the curvature of the unit-speed loops in (see [58] ) . To do this we use submanifold 
geometry. 

As in finite dimensional geometry, if A/" is a submanifold of a (weak) Riemannian 
manifold M then the induced Levi-Civita connection on N is related to that on M 
via the second fundamental form 

(5.1) Ii{U,V) = V^V -V^V 

where U and V are vector fields tangent to N . 11 is symmetric and tensorial (i.e., its 
value at any p G M depends only on the values Up and Vp) and the curvature of 
can then be computed using the Gauss-Codazzi formula 

{{R^{U, V)V, U)) = {{R^{U, V)V, U)) 

+ mu,U),Il{V,V))) - mU,V)MU,V))). 

li N = Diff^(M) is the volumorphism group and M = Diff(M) is the group of 
all diffeomorphisms with the metric then the corresponding second fundamental 
form is n(C/, V) = VA^^ div (Vnw) o r/ where U = uorj, V = voi] and where u and v 
are divergence free vector fields on M. The following theorem summarizes the known 
results in this important case; see [STJ |38l HSl |55l EOl Elj 

Theorem 5.1. Let M be a compact manifold of dimension n > 2, possibly with 
boundary. Consider the volumorphism group Diff^(M) with the metric. For any 
u G reDiff^(M) define 

( ^ ■ r {{R{v,u)u,v))l2 



min\'-^J — -li-^-^ , , II II 9 II II 9 // \\ 9 

i,eT,Diff^{M) ||w||2^2||it|li2 - {{u,v))l2 

to be the minimum sectional curvature in directions containing u, and similarly define 
Kmax{u) to be the maximum curvature. Then we have 

(1) iCmin(w) < 0, unless u is a Killing field in which case K^i^iiu) = 0; 

(2) if n >3 and M is flat then i^max(^^) > 0, unless div (V„m) = in which case 
i^max(li) = 0; 

(3) ifn = 2 and M is flat then i^max(^i) > 0, unless u is plane parallel u = f{x) 
or purely rotational u = f{r) M^, in either of which cases ^max(^^) = 0. 



Proof. All three statements follow from the Gauss-Codazzi formula (5.2) which is 
the most effective way to determine the sign of the curvature of Diff^(M). Part (1) 
can be found in [60j when n = 3 and M is flat although the technique works in any 
dimension. In part (2), the fact that div (V„ti) = implies i^max('u) = is essentially 
due to |45j while the converse can be proved using the same approximation scheme 
as in [56] (used there to find conjugate points). Finally, part (3) is a special case of 



a result in [55] which works for any steady flow on a surface. □ 
In particular, it follows from Theorem |5 . 1 1 that for any M of dimension n>2 the 



sectional curvature of Diff^(M) assumes both signs. The first examples in the 



special case of the flat 2-torus M = T were worked out by Arnold [Ij, who used the 



Lie-theoretic approach of Section 2.1 He derived a general formula for the sectional 



curvature of a group Ai = G with a right-invariant metric in terms of the coadjoint 
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operator (2.2) as 

{{R{u,v)v,u)) = -||ad*M + ad*w||^ — ((ad* u, ad*w)) 
' 3 , 1 



^,,ad„t;|p + ^((ad„i;,ad> - ad»^ 



and applied it to the case of Diff^ (T^ 



Example 5.2. If u and v are vector fields on with the stream functions f{x,y) = 
cos {jx + ky) and g{x, y) = cos {Ix + my) then the (unnormalized) sectional curvature 
of Difr^(T2) is 

, , „ , , „ TT'^{jm-kl)'^if + k'^ + P +m'^) 

(5.4) = - (o^,)^^^(,^i,;^)((^_,).^(J„,p) < 0. 

On the other hand, if we pick /(x, y) = cos (3A;x — y) + cos {2>kx + 2y) and g{x^ y) = 
cos {kx + y) + cos {kx — 2y) then 

9 

lim K{u,v) = 

Thus it is easy to find negative curvature, but there are many sections with positive 
curvature as well. 

There are two cases in which the curvature is known to have a remarkably simple 
form. The first is Diff(S'^) with the right-invariant metric whose curvature at the 
identity (and hence everywhere by right invar iance) is given by 

(5.5) m(u,.>..uU,^fiu.,-.u^f,. 

and thus is non- negative. We will show in Section [7] that this attractive formula does 
not generalize to higher dimensions and the sectional curvature of the right-invariant 
metric on Diff(M) can assume both signs. 

The second case is the space of densities Diff(M)/Diff^(M) equipped with the 



homogeneous Sobolev metric obtained by setting 6=1 and a = c = in (1.1). 
This space turns out to be isometric to the round sphere of radius 2 and therefore 
has constant positive curvature 

(5.6) {{R{u,v)v,u))^^ = 11^11^1 - {{u,v)f^,y 

We refer to |34| [29] for detailed calculations. 

When the exponential map is not smooth, the curvature may be positive but un- 
bounded above; this allows for conjugate points that occur arbitrarily close to any 
given point along a geodesic. In such situations one cannot determine stability study- 



ing geodesic deviation even for short times since the Rauch Theorem 4.2 immediately 
fails. This applies for example to the right-invariant metric on the Virasoro group 
(which yields the KdV equation) whose exponential map is also known not to be 
smooth and whose sectional curvature is unbounded and of both signs [461 llOj . 

Remark 5.3. In what follows we will often use the letter 5" as a shorthand notation 
for the (non-normalized) sectional curvature S{u,v) = {{R{u,v)v,u)) if the metric 
used is clear from the context. 
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6. The sign of the curvature: the one dimensional case 

In this section we present new results on the sign of the sectional curvature of the 
right-invariant metric on the group Diff(S'^). In this case the Sobolev H^-meticic 
(1.1) reduces to the a-b metric 

(6.1) {{u,v))fji = / [auv + buxVx) dx 



with a > and 6 > 0. The corresponding Euler- Arnold equation (2.1) reads 
(6.2) mt = —Sauux + b{2uxUxx + uuxxx), fn = Au = au — buxx- 

For a = 6 = 1 we get the periodic Camassa-Holm equation with period 1. For other 
values we can rescale by y = and s = t-^/f- so that (6.2) becomes 

Ug ^Syy ~t~ ^UUy '2UyUyy UUyyy 0, 

which is the Camassa-Holm equation with period \J ajb. 

Recall from Section [s] that the sectional curvature of Diff(5^) equipped with the a-b 
metric (6.1) is already known in the "end-point" cases where either a = or 6 = 0. 
In the former case, the sectional curvature is positive and constant (see (5.6)), and 
in the latter it is non- negative (see ( |5.5[ )). We will next show that when both a and 
b are positive, the curvature of Diff(S'^) can assume both signs. 

Lemma 6.1. The sectional curvature of Diff(S^) endowed with the right-invariant 
a-b metric (|6.1|) where a > and b > is given by 



(6.3) S{u,v) = {{R{u,v)v,u))hi = {{T{u,v),T{u,v)))h^ - {{T{u,u),T{v,v))) hi, 
for any u,v £ TeDiff(5^), where T is the Christoffel map defined by 



(6.4) T{u,v) = A ^dx ^auv -\- ^UxV^j^ , A = a — bdx- 

Proof. We have ad*u = A~^{2auvx + avux — 2bvxUxx — bvuxxx) and therefore an easily 
verified identity 

ad*ti + ad*^ = dx{uv) + 2T{u,v), 

which yields 

ad*ti = Idxiu^) + r(ii, u). 



Using the general curvature formula (5.3), we obtain 

{{R{u,v)v,u))hi = ||r(u,t;)||^i - {{r{u,u),T{v,v)))Hi + L{u,v), 

where 

L{u,v) = {{T{u,v),dx{uv)))H^ - l{{Tiu,u),dx{v^)))m - lmv,v),dx{u^)))m 

+ \\muv)\\j,^-l{{dx{n'),dx{v')))H^ 
- |||adtiu||^i + 2((aduU,ad*u - adlv))^! 

and aduV = —VxU-\-UxV. Since ((/, A~'^g))}ji = f^i fg dx for any functions / and g, we 
can perform all these computations without ever explicitly evaluating A^^ . A lengthy 
computation involving integration by parts shows that L(n, v) is always zero. □ 

The following proposition shows that it is easy to find sections of positive 
curvature on Diff(S'^) with the a-b metric. In fact, the curvature is strictly positive 
along all subspaces spanned by two trigonometric functions. 
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Theorem 6.2. Consider Diff(5^) endowed with the right-invariant metric given 
at the identity by (6.1) with a > and b > 0. If k and I are strictly positive distinct 
integer multiples of 2tt, then 

(6.5) S {cos kx , cos Ix) = S{cos kx,smlx) = S{smkx,sinlx) = C{k,l) > 0, 
where S{u,v) = {{R{u,v)v,u)) fji and 



C{k,l) 



1 / {a+^klf 2 {a-^klf 



;{k-ir + 



8 \a + b{k-iy' ' a + b{k + l)^ 
Moreover, for k an integer multiple of 2tt, we have 



{k + lf 



S{cos kx, sin kx) = 2C{k, k) 



a + Abk'^ 

S{cos kx, 1) = S{s\n kx, 1) = 2C{k, 0) = 



2u2 



a^k 



2(a + 6fc2) 



> 0. 



Proof. Let u = coskx and v = cos/x. From the definition (6.4) of the Cliristoffel map 
r we have 



r(u, u) 
T{u,v) 



sin {2kx), 



2ak - bk^ 
' 2a + 86A;2 
{k + l){2a - bkl) 
Aa + Ab{k + IY 



T{v,v) 



2al - bfi 
'2a + 8bP 



sin {2lx), 



. {k-l){2a + bkl) . 

sm ik + l)x — sm ik - l)x. 

Aa + Ab{k-IY 



Substituting these into (6.3) gives the formula for 5(cos /cx, cos /x). The other formu- 
las are proved in a similar way. □ 

Sections of negative curvature are trickier to find. 

Theorem 6.3. For every choice of a > and b > 0, there exist velocity fields u and 
V such that the sectional curvature ofDiS{S^) endowed with the a-b metric (6.1) is 
strictly negative, i.e., S{u,v) = {{R{u,v)v,u)) < 0. 



Proof. Set a = a/{A7r'^b) and pick 



U[X) 



where 



+ cos Attx, 



3 a' 



v{x) = sin 2ttx, 
a-2 



2 a{a + 4.) ' 



substitution into (6.3) yields 

S{u, v) 



267r^(a^ + IBa^ + 3570,2 _ 20q, _ 35) 



(Q + 9)(a + 4)2 

Observe that this quantity is negative for < a < 0.34. 

The second example is constructed differently and works when a > 0.34. Choose 
a positive integer j such that |y^a/0.34 < j < y^Q/0.34 and define r = a/j"^ so that 
0.34 < r < 1.36. Set 



(73r2 - 188r + 45)(r + 16) 



' 128(r + 9)(r - 2)2 

It is easy to see that this is defined in the range specified above. Set 

u(x) = cos 27rjx + -0 cos 47rjx and r;(x) = sin27rjx + 2'0sin47rjx. 
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Substituting into (6.3), we obtain 



S{u, v) 



P{r) 



64 (r + 9)2(r + 4)(r-2)2' 

where P{r) = 1435r6 + 21940r5 - 55074r^ - 222512r3 + 584323r2 - 215364r + 15552. 
This quantity S{u, v) is negative for 0.34 < r < 1.36, as desired. □ 

Remark 6.4 {The fiCH equation). An interesting example of a right-invariant Re- 
type metric on Diff(S'^) is given at the identity by 

{{u,v))m= cn{u)fj,{v) + / u {x)v (x) dx , 

for any positive constant c, where fJ-{u) := J^i u[x)dx is the mean value of the field 
over the circle. This metric yields yet another integrable evolution equation 



2cfJ,{u)Ux + 2UxUxx + UUx 







as a geodesic equation on the diffeomorphism group which "interpolates" between 
the Hunter-Saxton and Camassa-Holm equations. (This equation is sometimes called 
the /iHS or fiCB. equation.) The group Diff(5'^) equipped with the metric above 
admits sections of negative curvature; e.g., S{u,v) < whenever 

u{x) = h COS (4:7rkx), v{x) = sin {2Trkx) , 

c 

and k is any nonzero integer. 



7. The sign of the curvature: higher dimensions 

In order to simplify the formulas we present the results for the case when M is the 
flat torus T" = M^/Z". We recall the general formula for the Sobolev metric (1.1) in 
the form 



(7.1) 



{{u,v))ui=a / {u,v)dfi + b / dv!' ■ 6v^ dfi + c / {dv!' , dv^) dfi 

Jjn Jjn Jjn 



and observe from (5.6) that when a 
ture of Diff(T' 



c = then the corresponding sectional curva- 



is strictly positive and constant. 
In this section, we show that in the general a-b-c case the sectional curvature of 
assumes both signs. The case when at least two of the parameters 

the case b = i 



( [TTD on Diff(T 

a, b, c are nonzero is treated in Section |7.1 
7.2, and the case a 



on the subgroup Diff^^exCT ) is in Section 7.3 



is treated in Section 
Ah 2D 



examples discussed below generalize naturally to higher dimensions. 
7.1. The //^-metric on Diff(T"): the EPDiff equation. In the case when all 



parameters a, b and c of the metric (7.1) are strictly positive the Euler- Arnold 



equation (|2.1|) is a multidimensional generalization of the Camassa-Holm equation. 

we obtain 



In the special case where a 
the EPDiff equation 



1 and the manifold is a flat torui 



Theorem 7.1. If M is the flat torus T" and at least two of the parameters a, b, c are 
nonzero in (7.1), then the curvature o/Diff(T") takes on both signs. 
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In general, EPDiff involves the rough Laplacian V* V rather than the Hodge Laplacian dS + 5d; 
these operators differ by a Ricci curvature term due to the Bochner-Weitzenbock formula. If the 



manifold is Einstein the EPDiff metric is a special case of ( 7.1 1 
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Proof. The formula for the coadjoint operator for the a-b-c metric (7.1) is 
(7.2) ad> = ((div v)Au + d{Au, v) + LadAu) , 

where Av = av^ + bd5v^ + c6dv^; see j29j . 

If we pick u = f{x)-^ and v = g{x)-^ then ( |7.2[ ) gives 

ad> = (a - 6(9^)"^ (a(2g^/ + /^g) - b{2g^f^^ + gfxxx)^- 

The value of c is irrelevant in this case, since dn^ = dv^ = 0. Thus, formula (5.5) 
yields examples with positive curvature whenever a 7^ 0. If a = but b is nonzero, 
then positive curvature directions exist since this space is isometric to a sphere 



Finally, if a and b are both nonzero, then Proposition 6.3 yields examples of negative 
curvature. 

To finish the proof we need negative-curvature examples when 6 = with a and c 
both nonzero, and when a = while b and c are both nonzero. We will present them 
for the two-dimensional flat torus T^. 

Let u = f{x)§y andu = g{x)§^ so that Au = {af-cfxx)§y and = {ag-bg^x)^- 
If a 7^ 0, then using (7.2) we find 



ad> = (a - bdl) ^ (^aff^ - cf^fxx^ ^, 

ad> = (a - cdly'^ (a{g^f + gf^) - c{gxfxx + gf^xa 
ad*t> = 0, 

ad> = (a - 6(9^)-i i'iagga: - b{2g^g^^ + ggxxx)) ^• 



d_ 
9y' 



We also have aduV = gfx-g^- If f{x) = g{x) = sin/cx for some k which is an integer 
multiple of 27r, then it is easy to see that 

k{a + ck'^) . ^, n 

ad,,u = — ., , „, sm. Ikx 75^, 

2(a + 46/c2) 

_ k{a + ck'^) . ^, o 
ad,,n = , , „ sm2/ex 75^, 



_ 3fc(a + bA;2) . ^ 
2(a + 46A;2) '9^' 

ad„ti = ^ sin2/cx 

Note that these formulas are valid if a = as well, as long as 6 7^ and c 7^ 0. 
Indeed, in this case all vectors have to be projected to the orthogonal complement of 
the harmonic fieldsj^but since 

smkxdxdy= / / s'm2kx dx dy = 0, 



/o Jo Jo Jo 

we see that all components are already orthogonal to the harmonic fields. 
Substitution into (5.3) now gives the formula 

_ k'^{7a^ - Sa'^bk'^ + bGa'^ck'^ + Uack% + 76c^k'^a + IGOc^k'^b) 
^"^'^^ ~ 32{a + 4ck^){a + 4bk'^) " 

^^The metric (7.1 1 on Diff(T^) with a = is degenerate and only defined on the homogeneous 
space Diff(T^)/T^. Hence everything is only defined modulo harmonic fields on T^. 
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In particular, when a = then 

S{u,v) = -5cA;^/16 

and when 6 = we obtain 



S{u, v) 



32(a + 4c/c2) 



In either case the sectional curvature is negative for any k ^ 0. All these examples 
work on T" as well, if x and y denote the first two variables of the coordinate system. 

□ 

7.2. The L^-metric on DifF(T"): the Burgers equation. If 6 = c = then 
the formula (7.1) reduces to the inner product and the corresponding geodesic 
equation (2.1) is the multi-dimensional Burgers equation 

ut + VuU + div(u)n + ^V{u, u) = 0, 

also called the template matching equation; see |23J. However, in contrast with the 
one-dimensional case, the curvature of the right-invariant L^-metric assumes both 
signs when n > 2. For simplicity we only prove the result for n = 2. 

Proposition 7.2. The sectional curvature o/Diff(T^) equipped with the right-invariant 
metric (7.1) with b = c = is given by the formula (5.5) for any u = f{x) ^ and 

V = g{x) in which case S{u,v) > 0. On the other hand, if u = sin (27rx) and 

V = sin^ (27rx) then S{u,v) < 0. 



Proof. In this case the operator ad*M defined by (2.2) has the form 

ad*u = u div v + {i^du^)^ + V(n, v) . 

Thus, when u = f{x)§^ and v = g{x)-§^ we get ad> = [2gxf + gfx) ^ which is 
the same formula as in the one-dimensional case and the first part of the proposition 
follows. 

Furthermore, ii w = g{x)^ then we compute 

ad^n = 0, 
ad> = ^ffx-§^, 
ad*w^ = {fxg + fgx)§^, 



&dl,w = gg^-^. 



Combining these formulas with aduW = —fgx-§j^ in (5.3), we get 

S{u,w) = aj^ (^If^g"^ -2ffxggx^ dx. 

Taking f{x) = sin (27rx) and g{x) = sin^ {2itx), we find S{u,w) = — 157r^/16. □ 

A similar consideration in the general case can be summarized as the following 
statement. 

Theorem 7.3. The sectional curvature o/Diff(T") equipped with the right-invariant 
metric (i.e. the a-b-c metric (7.1) with a = I and b = c = 0) assumes both signs. 
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7.3. The homogeneous if^-metric on Difffj exCT^). Consider next the Lie group 
Diff^^exCT^) of exact volumorphisms of the flat torus which consists of symplectic 
diffeomorphisms preserving the center of mass. Its Lie algebra consists of Hamiltonian 
vector fields u = sgrad/ with / G C°°(T^). Following Arnold [I] we calculate the 
sectional curvature of the metric 

(7.3) {{u,v)) = c I {du\dv!')dn, 

on this group. It turns out to be just as convenient to work with a more general 
right-invariant metric given at the identity by 



(7.4) ((n,w)) = ((sgrad/, sgrad 5)) = / fAgd^i, 

where the (positive-definite, symmetric) operator defining the inner product is given 
by the formula A = A(A) for some function A: — >• M"*"!^ For a vector p G M^, we 
will write F{p) = A(|pp) for convenience. The metric ( |7.3[ ) corresponds to X{z) = cz^ 
and A = cA^. 

Theorem 7.4. Suppose f{x,y) = cos (jx + fey) and g{x,y) = cos {lx + my), where 
j, k, I, m are integer multiples of 2it. Set p = (j, k) and q = {I, m) and let u = sgrad / 
and V = sgrad g. Ther^ 



(7,5) S(u,.) = ^{K^W -F(,))'(^ + ^ 

-\[F{j, + q) + F(p - «)) + F(p) + F(q) 

where p Aq = jm — kl. 

Proof. RecaU that sgrad / = - fy§^ + fx§^ so that if {f,g] = fxgy- fygx denotes the 
Poisson bracket, then [sgrad/, sgrad 5] = sgrad{f, g}. Given any smooth functions /, 
g and h on the torus let u = sgrad f,v = sgrad 5 and w = sgrad h be the corresponding 
skew-gradients. Integration by parts gives 

{{adlu,w)) = - [ {Af){g,h}dfi 



/ h{g,Af}dfi = {{w,sgmdA ^{g,Af})), 
7t2 



from which we deduce that ad*n = sgrad A~^{g,Af}. 

We can assume that \p\, \q\, \p — q\, and + are all nonzero. Furthermore, since 
^/ = IpP/ s-iid ^9 = \q\'^9 fiiid that 
ad*ti = ad*t; = 0, 

ad*ti = —F(p) sgrad A^"'^6' and ad* w = F(g) sgrad A^^0, 

where 

G{x,y) = {f,g}{x,y) = ^{jm - kl)(^cos {{j - l)x + {k - m)y) 

— cos ((j -|- l)x + {k + m)y) 



More precisely, if a function / is written in an eigenbasis of the positive-definite Laplacian as 
/ = a,k(j)k with A(?!)fc = 7fe(?!)fc, then A/ = akX{'yk)(t>k- 



13a 



Note that for F{p) = \p\ formula (|7.5[) reproduces (5.4 1, up to a rescaling factor. 
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Clearly, we have = ^{p /\ q){^ — ip) where ip and •(/' are eigenfunctions satisfying 
Aip = F{p — q)ip and = F{p + q)ip. Combining the above formulas we obtain 

p Aq 



aduV 



ad*ii 



sgrad 9 
{p A q)F{p) 



2 

sgrad 



sgrad((^ — Ip), 

If Ip 



ad:.= (^^Aj^,grad 



F{p - q) F{p + , 

V i' 

F{p - q) F{p + q) 



Since, on the other hand, we have 

((sgrad (/?, sgrad (/?)) = ^F{p — q), ((sgrad 99, sgrad ^ 
and ((sgrad ip, sgrad ip)) = ^F{p + q), 



substituting into the sectional curvature formula (5.3) yields (7.5). 



□ 



Corollary 7.5. In the particular case when the metric is given by (7.3) the sectional 
curvature can assume both signs depending on p and q. 



Indeed, for the metric (7.3) we have F(p) = X(\p\ ) = c\p\ and a straightforward 
computation gives S{u, v) > when p = (lOvr, 0) and q = (Svr, 27r) while S{u, v) < 
when p = {2tt, 0) and q = (0, 27r). 

Appendix A. Semi-direct products and curvature formulas 



Formula (6.3) for the sectional curvature of the metric derived in Lemma 



6.1 



of Section [6] resembles the formula for the curvature of a Riemannian submanifold Af 
isometrically immersed in an ambient manifold A4, given by (5.2). Observe however 
that the sign of (5.2) is the opposite of the sign for (6.3). The next proposition shows 
how to rearrange (6.3) to make the analogy work. 



Proposition A.l. The sectional curvature (6.3) o/Diff(S'"'^) with the metric (6.1) 
can be rewritten in either of the two following forms 

{uv, - vu^f dx + lliT{u, u),T{v, v))) - ^{{T{u, v),T{u, v))) 

,2 



(A.l) S{u,v 
(A.2) S{u,v 



{a{uvx - vux) + ^{vxUxx - UxVxx)) dx 

+ ^((Q(n, u), Q{v, v))) - liiQiu, v), Q{u, v))), 
where the bilinear maps T and Q are defined by 

(A. 3) T{u,v) = A'^ {auv + \uxVx^ , Q{u,v) = —d'^A~^ {auv + ^UxVx) ■ 

Proof. Note that r(n, v) defined by ( |6.4[ ) is related to Q and T by r(n, v) = dxT{u, v) 
and Q{u,v) = —dxT{u,v). 

We prove ( |A.l ) first. Let u and v he any vector fields and set q = auv + ^UxVx- 
Then T{u,v) = A^^q, so that 



{{Tiu,v),Tiu,v))) -- 
Now using the identity A~^d'^ 
{{T{u,v),T{u,v))) = l 



dxT{u, v)AdxT{u, v) dx 
= -i + tA-\ we find 



/ qA ^d^qdx. 



[auv + ^UxVx) dx - -{{T{u,v),T{u,v))). 
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Using the same trick also on the other term in (6.3) gives (A.l). A similar technique 
reduces to dol). □ 



The formulas (A.l) and (A. 2) suggest that Diff(S'^) with the a-b metric (6.1) can 
be realized as an isometrically immersed submanifold of some simpler manifold. For 



(A. 2) such an immersion is described in the following theorem. 



Theorem A. 2. Let Diff(5^) >< C°^(5^) denote the semidirect product of the diffeo- 
morphism group Diff(S'^) with C^{S^) with group structure given by 

{v,F)-{^,G) = {vo^,Fo^ + G). 

Define a right-invariant Riemannian metric on this group by the inner product at 
the identity 



(A.4) 



{{{u,f),{v,g))) 



SI 



[auv + bfg)dx. 



Let T: Diff(S'i) ^ DiS{S^) t< C°°{S^) denote the map T(r/) = (??,lnr/^). Then T is 
an embedding and a group homomorphism and the right-invariant metric induced on 



Diff(S'-^) by (A.4) is the a-b metric (6.1). Furthermore, the curvature formula (A. 2) 



is the Gauss- Codazzi formula for the embedding T. 

Proof. The Lie algebra for G = Diff(S'^) x C°^(5'^) is the semidirect product q = 
Vect(S'^) >< C°°{S^), which has an interpretation as the space of first-order differential 
operators vD + / for v G Vect and / G C°°. For details on the geometry of the 
semidirect product see e.g., |391I67] . For this semi-direct product the adjoint operator 
has the form 

ad(„j)(t;,fi() = {-uv^ + UxV,vfx 



ugx 



which implies that the Arnold operator (|2.2|) is 
(A.5) 



2uxv + uvx + ^gfx, gux + gxu) . 



The general curvature formula (5.3) implies, after some simplifications, that the cur- 



vature of the semidirect product Diff(S'^) k C°°(S'^) can be written as 



(A.6) {{Riiu,f),iv,g))iv,g),iu,f))) 



[UV, 



+ i^iafx - f9x)f dx 



+ 4 / {gUx - fvx){gux - fvx + 8vfx - 8ugx) dx. 



Observe that T is a group homomorphism: 

T(r? o e) = (r? o ln{7jx o + H^x)) = T(r?) • T(e) 

using the chain rule. Smoothness of this map can be proved using Sobolev topol- 
ogy. For our purposes it is enough to note that T is formally an immersion since at 
the identity we have 

DTe{u) = {u,Ux), 

which is obviously injectivej^ Hence the induced metric on Diff(S'^) is given at the 
identity by 



{{u,v)) = {{DTiu),DTiv))) 



SI 



[auv + buxVx)dx, 



which is precisely the a-b metric (6.1). 



-'^^In the interpretation via differential operators, we consider operators of the form uD + u^ — Du. 
The homomorphism DT : uD Du is evidently a homomorphism of Lie algebras Vect — >■ Vect ix C°° . 
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The orthogonal complement of the image of -DTg in the metric (A.4) consists of 
those vectors of the form {^hx , h) for some function h: S*^ — )• M. We can now compute 
the second fundamental form (5.1 ) of the embedding: because of the right-invariance, 
we have {{Il{u,u),w)) = ((ad*u, tf)) for w G (Im DTe)^ where ad* is the operator 
given in ( A.5| ). We then have 

n(n, u) = {\hx, h) where h = -A~^dl {au^ + ju^) . 

Polarization yields Ii{u,v) = [^Q{u,v)x,Q{u,v)^ , where Q is given by (A. 3). Sub- 
stituting (A.6) when f = u' and g = v' , together with the second fundamental form, 
into (5.2) reproduces (A. 2). □ 



We conjecture that there is another embedding-homomorphism of Diff(S'^) into an 
infinite-dimensional Lie group with right-invariant metric which reproduces formula 



(A.l), but we do not know what it might be. 
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